home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / cblas / test_tpsv.c < prev    next >
Encoding:
C/C++ Source or Header  |  2002-04-18  |  43.4 KB  |  1,660 lines
  1. #include <gsl/gsl_test.h>
  2. #include <gsl/gsl_ieee_utils.h>
  3. #include <gsl/gsl_math.h>
  4. #include <gsl/gsl_cblas.h>
  5.  
  6. #include "tests.h"
  7.  
  8. void
  9. test_tpsv (void) {
  10. const double flteps = 1e-4, dbleps = 1e-6;
  11.   {
  12.    int order = 101;
  13.    int trans = 111;
  14.    int uplo = 121;
  15.    int diag = 131;
  16.    int N = 2;
  17.    float A[] = { -0.381f, 0.53f, 0.451f };
  18.    float X[] = { 0.144f, 0.032f };
  19.    int incX = -1;
  20.    float x_expected[] = { 0.31929f, 0.360168f };
  21.    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  22.    {
  23.      int i;
  24.      for (i = 0; i < 2; i++) {
  25.        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1310)");
  26.      }
  27.    };
  28.   };
  29.  
  30.  
  31.   {
  32.    int order = 101;
  33.    int trans = 111;
  34.    int uplo = 121;
  35.    int diag = 132;
  36.    int N = 2;
  37.    float A[] = { -0.381f, 0.53f, 0.451f };
  38.    float X[] = { 0.144f, 0.032f };
  39.    int incX = -1;
  40.    float x_expected[] = { 0.144f, -0.04432f };
  41.    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  42.    {
  43.      int i;
  44.      for (i = 0; i < 2; i++) {
  45.        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1311)");
  46.      }
  47.    };
  48.   };
  49.  
  50.  
  51.   {
  52.    int order = 101;
  53.    int trans = 111;
  54.    int uplo = 122;
  55.    int diag = 131;
  56.    int N = 2;
  57.    float A[] = { -0.381f, 0.53f, 0.451f };
  58.    float X[] = { 0.144f, 0.032f };
  59.    int incX = -1;
  60.    float x_expected[] = { 0.417992f, -0.0839895f };
  61.    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  62.    {
  63.      int i;
  64.      for (i = 0; i < 2; i++) {
  65.        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1312)");
  66.      }
  67.    };
  68.   };
  69.  
  70.  
  71.   {
  72.    int order = 101;
  73.    int trans = 111;
  74.    int uplo = 122;
  75.    int diag = 132;
  76.    int N = 2;
  77.    float A[] = { -0.381f, 0.53f, 0.451f };
  78.    float X[] = { 0.144f, 0.032f };
  79.    int incX = -1;
  80.    float x_expected[] = { 0.12704f, 0.032f };
  81.    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  82.    {
  83.      int i;
  84.      for (i = 0; i < 2; i++) {
  85.        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1313)");
  86.      }
  87.    };
  88.   };
  89.  
  90.  
  91.   {
  92.    int order = 102;
  93.    int trans = 111;
  94.    int uplo = 121;
  95.    int diag = 131;
  96.    int N = 2;
  97.    float A[] = { -0.381f, 0.53f, 0.451f };
  98.    float X[] = { 0.144f, 0.032f };
  99.    int incX = -1;
  100.    float x_expected[] = { 0.31929f, 0.360168f };
  101.    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  102.    {
  103.      int i;
  104.      for (i = 0; i < 2; i++) {
  105.        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1314)");
  106.      }
  107.    };
  108.   };
  109.  
  110.  
  111.   {
  112.    int order = 102;
  113.    int trans = 111;
  114.    int uplo = 121;
  115.    int diag = 132;
  116.    int N = 2;
  117.    float A[] = { -0.381f, 0.53f, 0.451f };
  118.    float X[] = { 0.144f, 0.032f };
  119.    int incX = -1;
  120.    float x_expected[] = { 0.144f, -0.04432f };
  121.    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  122.    {
  123.      int i;
  124.      for (i = 0; i < 2; i++) {
  125.        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1315)");
  126.      }
  127.    };
  128.   };
  129.  
  130.  
  131.   {
  132.    int order = 102;
  133.    int trans = 111;
  134.    int uplo = 122;
  135.    int diag = 131;
  136.    int N = 2;
  137.    float A[] = { -0.381f, 0.53f, 0.451f };
  138.    float X[] = { 0.144f, 0.032f };
  139.    int incX = -1;
  140.    float x_expected[] = { 0.417992f, -0.0839895f };
  141.    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  142.    {
  143.      int i;
  144.      for (i = 0; i < 2; i++) {
  145.        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1316)");
  146.      }
  147.    };
  148.   };
  149.  
  150.  
  151.   {
  152.    int order = 102;
  153.    int trans = 111;
  154.    int uplo = 122;
  155.    int diag = 132;
  156.    int N = 2;
  157.    float A[] = { -0.381f, 0.53f, 0.451f };
  158.    float X[] = { 0.144f, 0.032f };
  159.    int incX = -1;
  160.    float x_expected[] = { 0.12704f, 0.032f };
  161.    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  162.    {
  163.      int i;
  164.      for (i = 0; i < 2; i++) {
  165.        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1317)");
  166.      }
  167.    };
  168.   };
  169.  
  170.  
  171.   {
  172.    int order = 101;
  173.    int trans = 112;
  174.    int uplo = 121;
  175.    int diag = 131;
  176.    int N = 2;
  177.    float A[] = { -0.381f, 0.53f, 0.451f };
  178.    float X[] = { 0.144f, 0.032f };
  179.    int incX = -1;
  180.    float x_expected[] = { 0.417992f, -0.0839895f };
  181.    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  182.    {
  183.      int i;
  184.      for (i = 0; i < 2; i++) {
  185.        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1318)");
  186.      }
  187.    };
  188.   };
  189.  
  190.  
  191.   {
  192.    int order = 101;
  193.    int trans = 112;
  194.    int uplo = 121;
  195.    int diag = 132;
  196.    int N = 2;
  197.    float A[] = { -0.381f, 0.53f, 0.451f };
  198.    float X[] = { 0.144f, 0.032f };
  199.    int incX = -1;
  200.    float x_expected[] = { 0.12704f, 0.032f };
  201.    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  202.    {
  203.      int i;
  204.      for (i = 0; i < 2; i++) {
  205.        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1319)");
  206.      }
  207.    };
  208.   };
  209.  
  210.  
  211.   {
  212.    int order = 101;
  213.    int trans = 112;
  214.    int uplo = 122;
  215.    int diag = 131;
  216.    int N = 2;
  217.    float A[] = { -0.381f, 0.53f, 0.451f };
  218.    float X[] = { 0.144f, 0.032f };
  219.    int incX = -1;
  220.    float x_expected[] = { 0.31929f, 0.360168f };
  221.    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  222.    {
  223.      int i;
  224.      for (i = 0; i < 2; i++) {
  225.        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1320)");
  226.      }
  227.    };
  228.   };
  229.  
  230.  
  231.   {
  232.    int order = 101;
  233.    int trans = 112;
  234.    int uplo = 122;
  235.    int diag = 132;
  236.    int N = 2;
  237.    float A[] = { -0.381f, 0.53f, 0.451f };
  238.    float X[] = { 0.144f, 0.032f };
  239.    int incX = -1;
  240.    float x_expected[] = { 0.144f, -0.04432f };
  241.    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  242.    {
  243.      int i;
  244.      for (i = 0; i < 2; i++) {
  245.        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1321)");
  246.      }
  247.    };
  248.   };
  249.  
  250.  
  251.   {
  252.    int order = 102;
  253.    int trans = 112;
  254.    int uplo = 121;
  255.    int diag = 131;
  256.    int N = 2;
  257.    float A[] = { -0.381f, 0.53f, 0.451f };
  258.    float X[] = { 0.144f, 0.032f };
  259.    int incX = -1;
  260.    float x_expected[] = { 0.417992f, -0.0839895f };
  261.    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  262.    {
  263.      int i;
  264.      for (i = 0; i < 2; i++) {
  265.        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1322)");
  266.      }
  267.    };
  268.   };
  269.  
  270.  
  271.   {
  272.    int order = 102;
  273.    int trans = 112;
  274.    int uplo = 121;
  275.    int diag = 132;
  276.    int N = 2;
  277.    float A[] = { -0.381f, 0.53f, 0.451f };
  278.    float X[] = { 0.144f, 0.032f };
  279.    int incX = -1;
  280.    float x_expected[] = { 0.12704f, 0.032f };
  281.    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  282.    {
  283.      int i;
  284.      for (i = 0; i < 2; i++) {
  285.        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1323)");
  286.      }
  287.    };
  288.   };
  289.  
  290.  
  291.   {
  292.    int order = 102;
  293.    int trans = 112;
  294.    int uplo = 122;
  295.    int diag = 131;
  296.    int N = 2;
  297.    float A[] = { -0.381f, 0.53f, 0.451f };
  298.    float X[] = { 0.144f, 0.032f };
  299.    int incX = -1;
  300.    float x_expected[] = { 0.31929f, 0.360168f };
  301.    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  302.    {
  303.      int i;
  304.      for (i = 0; i < 2; i++) {
  305.        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1324)");
  306.      }
  307.    };
  308.   };
  309.  
  310.  
  311.   {
  312.    int order = 102;
  313.    int trans = 112;
  314.    int uplo = 122;
  315.    int diag = 132;
  316.    int N = 2;
  317.    float A[] = { -0.381f, 0.53f, 0.451f };
  318.    float X[] = { 0.144f, 0.032f };
  319.    int incX = -1;
  320.    float x_expected[] = { 0.144f, -0.04432f };
  321.    cblas_stpsv(order, uplo, trans, diag, N, A, X, incX);
  322.    {
  323.      int i;
  324.      for (i = 0; i < 2; i++) {
  325.        gsl_test_rel(X[i], x_expected[i], flteps, "stpsv(case 1325)");
  326.      }
  327.    };
  328.   };
  329.  
  330.  
  331.   {
  332.    int order = 101;
  333.    int trans = 111;
  334.    int uplo = 121;
  335.    int diag = 131;
  336.    int N = 2;
  337.    double A[] = { 0.355, -0.698, -0.42 };
  338.    double X[] = { -0.702, -0.661 };
  339.    int incX = -1;
  340.    double x_expected[] = { 1.67142857143, 1.42438631791 };
  341.    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
  342.    {
  343.      int i;
  344.      for (i = 0; i < 2; i++) {
  345.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1326)");
  346.      }
  347.    };
  348.   };
  349.  
  350.  
  351.   {
  352.    int order = 101;
  353.    int trans = 111;
  354.    int uplo = 121;
  355.    int diag = 132;
  356.    int N = 2;
  357.    double A[] = { 0.355, -0.698, -0.42 };
  358.    double X[] = { -0.702, -0.661 };
  359.    int incX = -1;
  360.    double x_expected[] = { -0.702, -1.150996 };
  361.    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
  362.    {
  363.      int i;
  364.      for (i = 0; i < 2; i++) {
  365.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1327)");
  366.      }
  367.    };
  368.   };
  369.  
  370.  
  371.   {
  372.    int order = 101;
  373.    int trans = 111;
  374.    int uplo = 122;
  375.    int diag = 131;
  376.    int N = 2;
  377.    double A[] = { 0.355, -0.698, -0.42 };
  378.    double X[] = { -0.702, -0.661 };
  379.    int incX = -1;
  380.    double x_expected[] = { 4.76584842388, -1.86197183099 };
  381.    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
  382.    {
  383.      int i;
  384.      for (i = 0; i < 2; i++) {
  385.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1328)");
  386.      }
  387.    };
  388.   };
  389.  
  390.  
  391.   {
  392.    int order = 101;
  393.    int trans = 111;
  394.    int uplo = 122;
  395.    int diag = 132;
  396.    int N = 2;
  397.    double A[] = { 0.355, -0.698, -0.42 };
  398.    double X[] = { -0.702, -0.661 };
  399.    int incX = -1;
  400.    double x_expected[] = { -1.163378, -0.661 };
  401.    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
  402.    {
  403.      int i;
  404.      for (i = 0; i < 2; i++) {
  405.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1329)");
  406.      }
  407.    };
  408.   };
  409.  
  410.  
  411.   {
  412.    int order = 102;
  413.    int trans = 111;
  414.    int uplo = 121;
  415.    int diag = 131;
  416.    int N = 2;
  417.    double A[] = { 0.355, -0.698, -0.42 };
  418.    double X[] = { -0.702, -0.661 };
  419.    int incX = -1;
  420.    double x_expected[] = { 1.67142857143, 1.42438631791 };
  421.    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
  422.    {
  423.      int i;
  424.      for (i = 0; i < 2; i++) {
  425.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1330)");
  426.      }
  427.    };
  428.   };
  429.  
  430.  
  431.   {
  432.    int order = 102;
  433.    int trans = 111;
  434.    int uplo = 121;
  435.    int diag = 132;
  436.    int N = 2;
  437.    double A[] = { 0.355, -0.698, -0.42 };
  438.    double X[] = { -0.702, -0.661 };
  439.    int incX = -1;
  440.    double x_expected[] = { -0.702, -1.150996 };
  441.    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
  442.    {
  443.      int i;
  444.      for (i = 0; i < 2; i++) {
  445.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1331)");
  446.      }
  447.    };
  448.   };
  449.  
  450.  
  451.   {
  452.    int order = 102;
  453.    int trans = 111;
  454.    int uplo = 122;
  455.    int diag = 131;
  456.    int N = 2;
  457.    double A[] = { 0.355, -0.698, -0.42 };
  458.    double X[] = { -0.702, -0.661 };
  459.    int incX = -1;
  460.    double x_expected[] = { 4.76584842388, -1.86197183099 };
  461.    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
  462.    {
  463.      int i;
  464.      for (i = 0; i < 2; i++) {
  465.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1332)");
  466.      }
  467.    };
  468.   };
  469.  
  470.  
  471.   {
  472.    int order = 102;
  473.    int trans = 111;
  474.    int uplo = 122;
  475.    int diag = 132;
  476.    int N = 2;
  477.    double A[] = { 0.355, -0.698, -0.42 };
  478.    double X[] = { -0.702, -0.661 };
  479.    int incX = -1;
  480.    double x_expected[] = { -1.163378, -0.661 };
  481.    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
  482.    {
  483.      int i;
  484.      for (i = 0; i < 2; i++) {
  485.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1333)");
  486.      }
  487.    };
  488.   };
  489.  
  490.  
  491.   {
  492.    int order = 101;
  493.    int trans = 112;
  494.    int uplo = 121;
  495.    int diag = 131;
  496.    int N = 2;
  497.    double A[] = { 0.355, -0.698, -0.42 };
  498.    double X[] = { -0.702, -0.661 };
  499.    int incX = -1;
  500.    double x_expected[] = { 4.76584842388, -1.86197183099 };
  501.    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
  502.    {
  503.      int i;
  504.      for (i = 0; i < 2; i++) {
  505.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1334)");
  506.      }
  507.    };
  508.   };
  509.  
  510.  
  511.   {
  512.    int order = 101;
  513.    int trans = 112;
  514.    int uplo = 121;
  515.    int diag = 132;
  516.    int N = 2;
  517.    double A[] = { 0.355, -0.698, -0.42 };
  518.    double X[] = { -0.702, -0.661 };
  519.    int incX = -1;
  520.    double x_expected[] = { -1.163378, -0.661 };
  521.    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
  522.    {
  523.      int i;
  524.      for (i = 0; i < 2; i++) {
  525.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1335)");
  526.      }
  527.    };
  528.   };
  529.  
  530.  
  531.   {
  532.    int order = 101;
  533.    int trans = 112;
  534.    int uplo = 122;
  535.    int diag = 131;
  536.    int N = 2;
  537.    double A[] = { 0.355, -0.698, -0.42 };
  538.    double X[] = { -0.702, -0.661 };
  539.    int incX = -1;
  540.    double x_expected[] = { 1.67142857143, 1.42438631791 };
  541.    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
  542.    {
  543.      int i;
  544.      for (i = 0; i < 2; i++) {
  545.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1336)");
  546.      }
  547.    };
  548.   };
  549.  
  550.  
  551.   {
  552.    int order = 101;
  553.    int trans = 112;
  554.    int uplo = 122;
  555.    int diag = 132;
  556.    int N = 2;
  557.    double A[] = { 0.355, -0.698, -0.42 };
  558.    double X[] = { -0.702, -0.661 };
  559.    int incX = -1;
  560.    double x_expected[] = { -0.702, -1.150996 };
  561.    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
  562.    {
  563.      int i;
  564.      for (i = 0; i < 2; i++) {
  565.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1337)");
  566.      }
  567.    };
  568.   };
  569.  
  570.  
  571.   {
  572.    int order = 102;
  573.    int trans = 112;
  574.    int uplo = 121;
  575.    int diag = 131;
  576.    int N = 2;
  577.    double A[] = { 0.355, -0.698, -0.42 };
  578.    double X[] = { -0.702, -0.661 };
  579.    int incX = -1;
  580.    double x_expected[] = { 4.76584842388, -1.86197183099 };
  581.    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
  582.    {
  583.      int i;
  584.      for (i = 0; i < 2; i++) {
  585.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1338)");
  586.      }
  587.    };
  588.   };
  589.  
  590.  
  591.   {
  592.    int order = 102;
  593.    int trans = 112;
  594.    int uplo = 121;
  595.    int diag = 132;
  596.    int N = 2;
  597.    double A[] = { 0.355, -0.698, -0.42 };
  598.    double X[] = { -0.702, -0.661 };
  599.    int incX = -1;
  600.    double x_expected[] = { -1.163378, -0.661 };
  601.    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
  602.    {
  603.      int i;
  604.      for (i = 0; i < 2; i++) {
  605.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1339)");
  606.      }
  607.    };
  608.   };
  609.  
  610.  
  611.   {
  612.    int order = 102;
  613.    int trans = 112;
  614.    int uplo = 122;
  615.    int diag = 131;
  616.    int N = 2;
  617.    double A[] = { 0.355, -0.698, -0.42 };
  618.    double X[] = { -0.702, -0.661 };
  619.    int incX = -1;
  620.    double x_expected[] = { 1.67142857143, 1.42438631791 };
  621.    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
  622.    {
  623.      int i;
  624.      for (i = 0; i < 2; i++) {
  625.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1340)");
  626.      }
  627.    };
  628.   };
  629.  
  630.  
  631.   {
  632.    int order = 102;
  633.    int trans = 112;
  634.    int uplo = 122;
  635.    int diag = 132;
  636.    int N = 2;
  637.    double A[] = { 0.355, -0.698, -0.42 };
  638.    double X[] = { -0.702, -0.661 };
  639.    int incX = -1;
  640.    double x_expected[] = { -0.702, -1.150996 };
  641.    cblas_dtpsv(order, uplo, trans, diag, N, A, X, incX);
  642.    {
  643.      int i;
  644.      for (i = 0; i < 2; i++) {
  645.        gsl_test_rel(X[i], x_expected[i], dbleps, "dtpsv(case 1341)");
  646.      }
  647.    };
  648.   };
  649.  
  650.  
  651.   {
  652.    int order = 101;
  653.    int trans = 111;
  654.    int uplo = 121;
  655.    int diag = 131;
  656.    int N = 2;
  657.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  658.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  659.    int incX = -1;
  660.    float x_expected[] = { -1.05533f, 0.0529057f, -3.93625f, 1.36003f };
  661.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  662.    {
  663.      int i;
  664.      for (i = 0; i < 2; i++) {
  665.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1342) real");
  666.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1342) imag");
  667.      };
  668.    };
  669.   };
  670.  
  671.  
  672.   {
  673.    int order = 101;
  674.    int trans = 111;
  675.    int uplo = 121;
  676.    int diag = 132;
  677.    int N = 2;
  678.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  679.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  680.    int incX = -1;
  681.    float x_expected[] = { -0.133f, 0.998f, 0.818576f, 0.163438f };
  682.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  683.    {
  684.      int i;
  685.      for (i = 0; i < 2; i++) {
  686.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1343) real");
  687.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1343) imag");
  688.      };
  689.    };
  690.   };
  691.  
  692.  
  693.   {
  694.    int order = 101;
  695.    int trans = 111;
  696.    int uplo = 122;
  697.    int diag = 131;
  698.    int N = 2;
  699.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  700.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  701.    int incX = -1;
  702.    float x_expected[] = { -2.48362f, 1.13085f, -1.67581f, -0.273264f };
  703.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  704.    {
  705.      int i;
  706.      for (i = 0; i < 2; i++) {
  707.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1344) real");
  708.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1344) imag");
  709.      };
  710.    };
  711.   };
  712.  
  713.  
  714.   {
  715.    int order = 101;
  716.    int trans = 111;
  717.    int uplo = 122;
  718.    int diag = 132;
  719.    int N = 2;
  720.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  721.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  722.    int incX = -1;
  723.    float x_expected[] = { 0.431924f, 0.679112f, -0.072f, 0.642f };
  724.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  725.    {
  726.      int i;
  727.      for (i = 0; i < 2; i++) {
  728.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1345) real");
  729.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1345) imag");
  730.      };
  731.    };
  732.   };
  733.  
  734.  
  735.   {
  736.    int order = 102;
  737.    int trans = 111;
  738.    int uplo = 121;
  739.    int diag = 131;
  740.    int N = 2;
  741.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  742.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  743.    int incX = -1;
  744.    float x_expected[] = { -1.05533f, 0.0529057f, -3.93625f, 1.36003f };
  745.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  746.    {
  747.      int i;
  748.      for (i = 0; i < 2; i++) {
  749.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1346) real");
  750.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1346) imag");
  751.      };
  752.    };
  753.   };
  754.  
  755.  
  756.   {
  757.    int order = 102;
  758.    int trans = 111;
  759.    int uplo = 121;
  760.    int diag = 132;
  761.    int N = 2;
  762.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  763.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  764.    int incX = -1;
  765.    float x_expected[] = { -0.133f, 0.998f, 0.818576f, 0.163438f };
  766.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  767.    {
  768.      int i;
  769.      for (i = 0; i < 2; i++) {
  770.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1347) real");
  771.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1347) imag");
  772.      };
  773.    };
  774.   };
  775.  
  776.  
  777.   {
  778.    int order = 102;
  779.    int trans = 111;
  780.    int uplo = 122;
  781.    int diag = 131;
  782.    int N = 2;
  783.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  784.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  785.    int incX = -1;
  786.    float x_expected[] = { -2.48362f, 1.13085f, -1.67581f, -0.273264f };
  787.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  788.    {
  789.      int i;
  790.      for (i = 0; i < 2; i++) {
  791.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1348) real");
  792.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1348) imag");
  793.      };
  794.    };
  795.   };
  796.  
  797.  
  798.   {
  799.    int order = 102;
  800.    int trans = 111;
  801.    int uplo = 122;
  802.    int diag = 132;
  803.    int N = 2;
  804.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  805.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  806.    int incX = -1;
  807.    float x_expected[] = { 0.431924f, 0.679112f, -0.072f, 0.642f };
  808.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  809.    {
  810.      int i;
  811.      for (i = 0; i < 2; i++) {
  812.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1349) real");
  813.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1349) imag");
  814.      };
  815.    };
  816.   };
  817.  
  818.  
  819.   {
  820.    int order = 101;
  821.    int trans = 112;
  822.    int uplo = 121;
  823.    int diag = 131;
  824.    int N = 2;
  825.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  826.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  827.    int incX = -1;
  828.    float x_expected[] = { -2.48362f, 1.13085f, -1.67581f, -0.273264f };
  829.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  830.    {
  831.      int i;
  832.      for (i = 0; i < 2; i++) {
  833.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1350) real");
  834.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1350) imag");
  835.      };
  836.    };
  837.   };
  838.  
  839.  
  840.   {
  841.    int order = 101;
  842.    int trans = 112;
  843.    int uplo = 121;
  844.    int diag = 132;
  845.    int N = 2;
  846.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  847.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  848.    int incX = -1;
  849.    float x_expected[] = { 0.431924f, 0.679112f, -0.072f, 0.642f };
  850.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  851.    {
  852.      int i;
  853.      for (i = 0; i < 2; i++) {
  854.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1351) real");
  855.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1351) imag");
  856.      };
  857.    };
  858.   };
  859.  
  860.  
  861.   {
  862.    int order = 101;
  863.    int trans = 112;
  864.    int uplo = 122;
  865.    int diag = 131;
  866.    int N = 2;
  867.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  868.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  869.    int incX = -1;
  870.    float x_expected[] = { -1.05533f, 0.0529057f, -3.93625f, 1.36003f };
  871.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  872.    {
  873.      int i;
  874.      for (i = 0; i < 2; i++) {
  875.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1352) real");
  876.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1352) imag");
  877.      };
  878.    };
  879.   };
  880.  
  881.  
  882.   {
  883.    int order = 101;
  884.    int trans = 112;
  885.    int uplo = 122;
  886.    int diag = 132;
  887.    int N = 2;
  888.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  889.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  890.    int incX = -1;
  891.    float x_expected[] = { -0.133f, 0.998f, 0.818576f, 0.163438f };
  892.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  893.    {
  894.      int i;
  895.      for (i = 0; i < 2; i++) {
  896.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1353) real");
  897.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1353) imag");
  898.      };
  899.    };
  900.   };
  901.  
  902.  
  903.   {
  904.    int order = 102;
  905.    int trans = 112;
  906.    int uplo = 121;
  907.    int diag = 131;
  908.    int N = 2;
  909.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  910.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  911.    int incX = -1;
  912.    float x_expected[] = { -2.48362f, 1.13085f, -1.67581f, -0.273264f };
  913.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  914.    {
  915.      int i;
  916.      for (i = 0; i < 2; i++) {
  917.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1354) real");
  918.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1354) imag");
  919.      };
  920.    };
  921.   };
  922.  
  923.  
  924.   {
  925.    int order = 102;
  926.    int trans = 112;
  927.    int uplo = 121;
  928.    int diag = 132;
  929.    int N = 2;
  930.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  931.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  932.    int incX = -1;
  933.    float x_expected[] = { 0.431924f, 0.679112f, -0.072f, 0.642f };
  934.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  935.    {
  936.      int i;
  937.      for (i = 0; i < 2; i++) {
  938.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1355) real");
  939.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1355) imag");
  940.      };
  941.    };
  942.   };
  943.  
  944.  
  945.   {
  946.    int order = 102;
  947.    int trans = 112;
  948.    int uplo = 122;
  949.    int diag = 131;
  950.    int N = 2;
  951.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  952.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  953.    int incX = -1;
  954.    float x_expected[] = { -1.05533f, 0.0529057f, -3.93625f, 1.36003f };
  955.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  956.    {
  957.      int i;
  958.      for (i = 0; i < 2; i++) {
  959.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1356) real");
  960.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1356) imag");
  961.      };
  962.    };
  963.   };
  964.  
  965.  
  966.   {
  967.    int order = 102;
  968.    int trans = 112;
  969.    int uplo = 122;
  970.    int diag = 132;
  971.    int N = 2;
  972.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  973.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  974.    int incX = -1;
  975.    float x_expected[] = { -0.133f, 0.998f, 0.818576f, 0.163438f };
  976.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  977.    {
  978.      int i;
  979.      for (i = 0; i < 2; i++) {
  980.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1357) real");
  981.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1357) imag");
  982.      };
  983.    };
  984.   };
  985.  
  986.  
  987.   {
  988.    int order = 101;
  989.    int trans = 113;
  990.    int uplo = 121;
  991.    int diag = 131;
  992.    int N = 2;
  993.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  994.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  995.    int incX = -1;
  996.    float x_expected[] = { 2.15867f, 1.69498f, 1.69471f, 0.104738f };
  997.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  998.    {
  999.      int i;
  1000.      for (i = 0; i < 2; i++) {
  1001.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1358) real");
  1002.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1358) imag");
  1003.      };
  1004.    };
  1005.   };
  1006.  
  1007.  
  1008.   {
  1009.    int order = 101;
  1010.    int trans = 113;
  1011.    int uplo = 121;
  1012.    int diag = 132;
  1013.    int N = 2;
  1014.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  1015.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  1016.    int incX = -1;
  1017.    float x_expected[] = { -0.613252f, 0.561896f, -0.072f, 0.642f };
  1018.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  1019.    {
  1020.      int i;
  1021.      for (i = 0; i < 2; i++) {
  1022.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1359) real");
  1023.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1359) imag");
  1024.      };
  1025.    };
  1026.   };
  1027.  
  1028.  
  1029.   {
  1030.    int order = 101;
  1031.    int trans = 113;
  1032.    int uplo = 122;
  1033.    int diag = 131;
  1034.    int N = 2;
  1035.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  1036.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  1037.    int incX = -1;
  1038.    float x_expected[] = { 1.00465f, 0.327432f, 3.44853f, 2.273f };
  1039.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  1040.    {
  1041.      int i;
  1042.      for (i = 0; i < 2; i++) {
  1043.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1360) real");
  1044.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1360) imag");
  1045.      };
  1046.    };
  1047.   };
  1048.  
  1049.  
  1050.   {
  1051.    int order = 101;
  1052.    int trans = 113;
  1053.    int uplo = 122;
  1054.    int diag = 132;
  1055.    int N = 2;
  1056.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  1057.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  1058.    int incX = -1;
  1059.    float x_expected[] = { -0.133f, 0.998f, -0.806168f, -0.053086f };
  1060.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  1061.    {
  1062.      int i;
  1063.      for (i = 0; i < 2; i++) {
  1064.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1361) real");
  1065.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1361) imag");
  1066.      };
  1067.    };
  1068.   };
  1069.  
  1070.  
  1071.   {
  1072.    int order = 102;
  1073.    int trans = 113;
  1074.    int uplo = 121;
  1075.    int diag = 131;
  1076.    int N = 2;
  1077.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  1078.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  1079.    int incX = -1;
  1080.    float x_expected[] = { 2.15867f, 1.69498f, 1.69471f, 0.104738f };
  1081.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  1082.    {
  1083.      int i;
  1084.      for (i = 0; i < 2; i++) {
  1085.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1362) real");
  1086.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1362) imag");
  1087.      };
  1088.    };
  1089.   };
  1090.  
  1091.  
  1092.   {
  1093.    int order = 102;
  1094.    int trans = 113;
  1095.    int uplo = 121;
  1096.    int diag = 132;
  1097.    int N = 2;
  1098.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  1099.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  1100.    int incX = -1;
  1101.    float x_expected[] = { -0.613252f, 0.561896f, -0.072f, 0.642f };
  1102.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  1103.    {
  1104.      int i;
  1105.      for (i = 0; i < 2; i++) {
  1106.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1363) real");
  1107.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1363) imag");
  1108.      };
  1109.    };
  1110.   };
  1111.  
  1112.  
  1113.   {
  1114.    int order = 102;
  1115.    int trans = 113;
  1116.    int uplo = 122;
  1117.    int diag = 131;
  1118.    int N = 2;
  1119.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  1120.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  1121.    int incX = -1;
  1122.    float x_expected[] = { 1.00465f, 0.327432f, 3.44853f, 2.273f };
  1123.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  1124.    {
  1125.      int i;
  1126.      for (i = 0; i < 2; i++) {
  1127.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1364) real");
  1128.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1364) imag");
  1129.      };
  1130.    };
  1131.   };
  1132.  
  1133.  
  1134.   {
  1135.    int order = 102;
  1136.    int trans = 113;
  1137.    int uplo = 122;
  1138.    int diag = 132;
  1139.    int N = 2;
  1140.    float A[] = { -0.019f, -0.38f, 0.588f, 0.814f, 0.173f, -0.937f };
  1141.    float X[] = { -0.133f, 0.998f, -0.072f, 0.642f };
  1142.    int incX = -1;
  1143.    float x_expected[] = { -0.133f, 0.998f, -0.806168f, -0.053086f };
  1144.    cblas_ctpsv(order, uplo, trans, diag, N, A, X, incX);
  1145.    {
  1146.      int i;
  1147.      for (i = 0; i < 2; i++) {
  1148.        gsl_test_rel(X[2*i], x_expected[2*i], flteps, "ctpsv(case 1365) real");
  1149.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], flteps, "ctpsv(case 1365) imag");
  1150.      };
  1151.    };
  1152.   };
  1153.  
  1154.  
  1155.   {
  1156.    int order = 101;
  1157.    int trans = 111;
  1158.    int uplo = 121;
  1159.    int diag = 131;
  1160.    int N = 2;
  1161.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1162.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1163.    int incX = -1;
  1164.    double x_expected[] = { 0.430509772467, -0.0927067365535, -0.611144484555, 0.999982608216 };
  1165.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1166.    {
  1167.      int i;
  1168.      for (i = 0; i < 2; i++) {
  1169.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1366) real");
  1170.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1366) imag");
  1171.      };
  1172.    };
  1173.   };
  1174.  
  1175.  
  1176.   {
  1177.    int order = 101;
  1178.    int trans = 111;
  1179.    int uplo = 121;
  1180.    int diag = 132;
  1181.    int N = 2;
  1182.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1183.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1184.    int incX = -1;
  1185.    double x_expected[] = { 0.344, -0.143, -0.795928, -0.523879 };
  1186.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1187.    {
  1188.      int i;
  1189.      for (i = 0; i < 2; i++) {
  1190.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1367) real");
  1191.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1367) imag");
  1192.      };
  1193.    };
  1194.   };
  1195.  
  1196.  
  1197.   {
  1198.    int order = 101;
  1199.    int trans = 111;
  1200.    int uplo = 122;
  1201.    int diag = 131;
  1202.    int N = 2;
  1203.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1204.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1205.    int incX = -1;
  1206.    double x_expected[] = { 1.00136463678, -1.84591534629, -1.12140892769, 0.696784840869 };
  1207.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1208.    {
  1209.      int i;
  1210.      for (i = 0; i < 2; i++) {
  1211.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1368) real");
  1212.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1368) imag");
  1213.      };
  1214.    };
  1215.   };
  1216.  
  1217.  
  1218.   {
  1219.    int order = 101;
  1220.    int trans = 111;
  1221.    int uplo = 122;
  1222.    int diag = 132;
  1223.    int N = 2;
  1224.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1225.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1226.    int incX = -1;
  1227.    double x_expected[] = { 1.707508, -0.042521, -0.668, -0.945 };
  1228.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1229.    {
  1230.      int i;
  1231.      for (i = 0; i < 2; i++) {
  1232.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1369) real");
  1233.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1369) imag");
  1234.      };
  1235.    };
  1236.   };
  1237.  
  1238.  
  1239.   {
  1240.    int order = 102;
  1241.    int trans = 111;
  1242.    int uplo = 121;
  1243.    int diag = 131;
  1244.    int N = 2;
  1245.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1246.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1247.    int incX = -1;
  1248.    double x_expected[] = { 0.430509772467, -0.0927067365535, -0.611144484555, 0.999982608216 };
  1249.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1250.    {
  1251.      int i;
  1252.      for (i = 0; i < 2; i++) {
  1253.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1370) real");
  1254.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1370) imag");
  1255.      };
  1256.    };
  1257.   };
  1258.  
  1259.  
  1260.   {
  1261.    int order = 102;
  1262.    int trans = 111;
  1263.    int uplo = 121;
  1264.    int diag = 132;
  1265.    int N = 2;
  1266.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1267.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1268.    int incX = -1;
  1269.    double x_expected[] = { 0.344, -0.143, -0.795928, -0.523879 };
  1270.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1271.    {
  1272.      int i;
  1273.      for (i = 0; i < 2; i++) {
  1274.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1371) real");
  1275.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1371) imag");
  1276.      };
  1277.    };
  1278.   };
  1279.  
  1280.  
  1281.   {
  1282.    int order = 102;
  1283.    int trans = 111;
  1284.    int uplo = 122;
  1285.    int diag = 131;
  1286.    int N = 2;
  1287.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1288.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1289.    int incX = -1;
  1290.    double x_expected[] = { 1.00136463678, -1.84591534629, -1.12140892769, 0.696784840869 };
  1291.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1292.    {
  1293.      int i;
  1294.      for (i = 0; i < 2; i++) {
  1295.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1372) real");
  1296.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1372) imag");
  1297.      };
  1298.    };
  1299.   };
  1300.  
  1301.  
  1302.   {
  1303.    int order = 102;
  1304.    int trans = 111;
  1305.    int uplo = 122;
  1306.    int diag = 132;
  1307.    int N = 2;
  1308.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1309.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1310.    int incX = -1;
  1311.    double x_expected[] = { 1.707508, -0.042521, -0.668, -0.945 };
  1312.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1313.    {
  1314.      int i;
  1315.      for (i = 0; i < 2; i++) {
  1316.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1373) real");
  1317.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1373) imag");
  1318.      };
  1319.    };
  1320.   };
  1321.  
  1322.  
  1323.   {
  1324.    int order = 101;
  1325.    int trans = 112;
  1326.    int uplo = 121;
  1327.    int diag = 131;
  1328.    int N = 2;
  1329.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1330.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1331.    int incX = -1;
  1332.    double x_expected[] = { 1.00136463678, -1.84591534629, -1.12140892769, 0.696784840869 };
  1333.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1334.    {
  1335.      int i;
  1336.      for (i = 0; i < 2; i++) {
  1337.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1374) real");
  1338.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1374) imag");
  1339.      };
  1340.    };
  1341.   };
  1342.  
  1343.  
  1344.   {
  1345.    int order = 101;
  1346.    int trans = 112;
  1347.    int uplo = 121;
  1348.    int diag = 132;
  1349.    int N = 2;
  1350.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1351.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1352.    int incX = -1;
  1353.    double x_expected[] = { 1.707508, -0.042521, -0.668, -0.945 };
  1354.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1355.    {
  1356.      int i;
  1357.      for (i = 0; i < 2; i++) {
  1358.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1375) real");
  1359.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1375) imag");
  1360.      };
  1361.    };
  1362.   };
  1363.  
  1364.  
  1365.   {
  1366.    int order = 101;
  1367.    int trans = 112;
  1368.    int uplo = 122;
  1369.    int diag = 131;
  1370.    int N = 2;
  1371.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1372.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1373.    int incX = -1;
  1374.    double x_expected[] = { 0.430509772467, -0.0927067365535, -0.611144484555, 0.999982608216 };
  1375.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1376.    {
  1377.      int i;
  1378.      for (i = 0; i < 2; i++) {
  1379.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1376) real");
  1380.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1376) imag");
  1381.      };
  1382.    };
  1383.   };
  1384.  
  1385.  
  1386.   {
  1387.    int order = 101;
  1388.    int trans = 112;
  1389.    int uplo = 122;
  1390.    int diag = 132;
  1391.    int N = 2;
  1392.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1393.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1394.    int incX = -1;
  1395.    double x_expected[] = { 0.344, -0.143, -0.795928, -0.523879 };
  1396.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1397.    {
  1398.      int i;
  1399.      for (i = 0; i < 2; i++) {
  1400.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1377) real");
  1401.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1377) imag");
  1402.      };
  1403.    };
  1404.   };
  1405.  
  1406.  
  1407.   {
  1408.    int order = 102;
  1409.    int trans = 112;
  1410.    int uplo = 121;
  1411.    int diag = 131;
  1412.    int N = 2;
  1413.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1414.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1415.    int incX = -1;
  1416.    double x_expected[] = { 1.00136463678, -1.84591534629, -1.12140892769, 0.696784840869 };
  1417.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1418.    {
  1419.      int i;
  1420.      for (i = 0; i < 2; i++) {
  1421.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1378) real");
  1422.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1378) imag");
  1423.      };
  1424.    };
  1425.   };
  1426.  
  1427.  
  1428.   {
  1429.    int order = 102;
  1430.    int trans = 112;
  1431.    int uplo = 121;
  1432.    int diag = 132;
  1433.    int N = 2;
  1434.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1435.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1436.    int incX = -1;
  1437.    double x_expected[] = { 1.707508, -0.042521, -0.668, -0.945 };
  1438.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1439.    {
  1440.      int i;
  1441.      for (i = 0; i < 2; i++) {
  1442.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1379) real");
  1443.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1379) imag");
  1444.      };
  1445.    };
  1446.   };
  1447.  
  1448.  
  1449.   {
  1450.    int order = 102;
  1451.    int trans = 112;
  1452.    int uplo = 122;
  1453.    int diag = 131;
  1454.    int N = 2;
  1455.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1456.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1457.    int incX = -1;
  1458.    double x_expected[] = { 0.430509772467, -0.0927067365535, -0.611144484555, 0.999982608216 };
  1459.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1460.    {
  1461.      int i;
  1462.      for (i = 0; i < 2; i++) {
  1463.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1380) real");
  1464.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1380) imag");
  1465.      };
  1466.    };
  1467.   };
  1468.  
  1469.  
  1470.   {
  1471.    int order = 102;
  1472.    int trans = 112;
  1473.    int uplo = 122;
  1474.    int diag = 132;
  1475.    int N = 2;
  1476.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1477.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1478.    int incX = -1;
  1479.    double x_expected[] = { 0.344, -0.143, -0.795928, -0.523879 };
  1480.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1481.    {
  1482.      int i;
  1483.      for (i = 0; i < 2; i++) {
  1484.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1381) real");
  1485.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1381) imag");
  1486.      };
  1487.    };
  1488.   };
  1489.  
  1490.  
  1491.   {
  1492.    int order = 101;
  1493.    int trans = 113;
  1494.    int uplo = 121;
  1495.    int diag = 131;
  1496.    int N = 2;
  1497.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1498.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1499.    int incX = -1;
  1500.    double x_expected[] = { -1.47384781823, -0.286556198408, 1.03098932879, -0.824698794397 };
  1501.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1502.    {
  1503.      int i;
  1504.      for (i = 0; i < 2; i++) {
  1505.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1382) real");
  1506.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1382) imag");
  1507.      };
  1508.    };
  1509.   };
  1510.  
  1511.  
  1512.   {
  1513.    int order = 101;
  1514.    int trans = 113;
  1515.    int uplo = 121;
  1516.    int diag = 132;
  1517.    int N = 2;
  1518.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1519.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1520.    int incX = -1;
  1521.    double x_expected[] = { -0.016172, 1.175911, -0.668, -0.945 };
  1522.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1523.    {
  1524.      int i;
  1525.      for (i = 0; i < 2; i++) {
  1526.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1383) real");
  1527.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1383) imag");
  1528.      };
  1529.    };
  1530.   };
  1531.  
  1532.  
  1533.   {
  1534.    int order = 101;
  1535.    int trans = 113;
  1536.    int uplo = 122;
  1537.    int diag = 131;
  1538.    int N = 2;
  1539.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1540.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1541.    int incX = -1;
  1542.    double x_expected[] = { 0.369363905801, -0.239798891331, 1.1759505739, -1.40027235656 };
  1543.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1544.    {
  1545.      int i;
  1546.      for (i = 0; i < 2; i++) {
  1547.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1384) real");
  1548.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1384) imag");
  1549.      };
  1550.    };
  1551.   };
  1552.  
  1553.  
  1554.   {
  1555.    int order = 101;
  1556.    int trans = 113;
  1557.    int uplo = 122;
  1558.    int diag = 132;
  1559.    int N = 2;
  1560.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1561.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1562.    int incX = -1;
  1563.    double x_expected[] = { 0.344, -0.143, -1.05676, -1.151335 };
  1564.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1565.    {
  1566.      int i;
  1567.      for (i = 0; i < 2; i++) {
  1568.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1385) real");
  1569.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1385) imag");
  1570.      };
  1571.    };
  1572.   };
  1573.  
  1574.  
  1575.   {
  1576.    int order = 102;
  1577.    int trans = 113;
  1578.    int uplo = 121;
  1579.    int diag = 131;
  1580.    int N = 2;
  1581.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1582.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1583.    int incX = -1;
  1584.    double x_expected[] = { -1.47384781823, -0.286556198408, 1.03098932879, -0.824698794397 };
  1585.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1586.    {
  1587.      int i;
  1588.      for (i = 0; i < 2; i++) {
  1589.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1386) real");
  1590.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1386) imag");
  1591.      };
  1592.    };
  1593.   };
  1594.  
  1595.  
  1596.   {
  1597.    int order = 102;
  1598.    int trans = 113;
  1599.    int uplo = 121;
  1600.    int diag = 132;
  1601.    int N = 2;
  1602.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1603.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1604.    int incX = -1;
  1605.    double x_expected[] = { -0.016172, 1.175911, -0.668, -0.945 };
  1606.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1607.    {
  1608.      int i;
  1609.      for (i = 0; i < 2; i++) {
  1610.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1387) real");
  1611.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1387) imag");
  1612.      };
  1613.    };
  1614.   };
  1615.  
  1616.  
  1617.   {
  1618.    int order = 102;
  1619.    int trans = 113;
  1620.    int uplo = 122;
  1621.    int diag = 131;
  1622.    int N = 2;
  1623.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1624.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1625.    int incX = -1;
  1626.    double x_expected[] = { 0.369363905801, -0.239798891331, 1.1759505739, -1.40027235656 };
  1627.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1628.    {
  1629.      int i;
  1630.      for (i = 0; i < 2; i++) {
  1631.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1388) real");
  1632.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1388) imag");
  1633.      };
  1634.    };
  1635.   };
  1636.  
  1637.  
  1638.   {
  1639.    int order = 102;
  1640.    int trans = 113;
  1641.    int uplo = 122;
  1642.    int diag = 132;
  1643.    int N = 2;
  1644.    double A[] = { 0.052, 0.875, 0.751, -0.912, 0.832, -0.153 };
  1645.    double X[] = { 0.344, -0.143, -0.668, -0.945 };
  1646.    int incX = -1;
  1647.    double x_expected[] = { 0.344, -0.143, -1.05676, -1.151335 };
  1648.    cblas_ztpsv(order, uplo, trans, diag, N, A, X, incX);
  1649.    {
  1650.      int i;
  1651.      for (i = 0; i < 2; i++) {
  1652.        gsl_test_rel(X[2*i], x_expected[2*i], dbleps, "ztpsv(case 1389) real");
  1653.        gsl_test_rel(X[2*i+1], x_expected[2*i+1], dbleps, "ztpsv(case 1389) imag");
  1654.      };
  1655.    };
  1656.   };
  1657.  
  1658.  
  1659. }
  1660.